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Abstract. We have constructed self-consistent light and mass distribution models for four disk galaxies at redshifts 
2 = 0.48, 0.58, 0.81 and 0.88, using the HST archive WFPC2 observations (HDF-S, MDS, Groth Strip survey) 
and rotation curves measured by Vogt et al. 119961 and Rigopoulou et al. (2002;). The models consist of three 
components: a bulge, a disk and a dark matter halo. Similarly to the sample studied in Paper I (Tamm & Tenjes, 
I2003II . light distribution of the galaxies in the outer parts is clearly steeper than a simple exponential disk. After 
applying k-corrections, calculated mass-to-light ratios for galactic disks within the maximum disk assumption are 
M/ Lb = 0.9, 7.4, 4.3 and 1.4, respectively. Together with the galaxies from Paper I, the mean {M/ Lb) ~ 2.5 at 
(z) ~ 0.9, indicating no significant evolution of M/ Lb with redshift. Central densities of dark matter halos for an 
isothermal model are 0.008, 0.035, 0.013, and 0.022 in units Mq/pc'^, respectively. Together with the galaxies from 
Paper I, the DM central density of the four galaxies at mean readshift (2:) ~ 0.9 is p(0) = (0.012 — 0.028) Mq/pc^, 
also showing no significant evolution with redshift. 

We have also constructed mass distribution models without assuming fiat rotation curves and without a dark 
matter component. Due to the limited extent of the observed rotation, the models without a dark halo nearly fit 
the observations. In this case, mass-to-light ratios for the galactic disks are 1.8, 9.7, 12. and 1.9, respectively. 

Key words, galaxies: photometry - galaxies: fundamental parameters - galaxies: high redshift ~ galaxies: spiral - 
galaxies: structure - dark matter 



ray cluster at z = 1.2 and in addition to the brightening 
they found also flattening of the faint end slope of the LF. 



^ 1. Introduction redshifts z ^ 5, PoH et al. H2UU3|I in the HDS-S up to z 

^, 1 n 1 1 11 1 1 • • 1 3.5, and certain evolutionary changes have been found. 

• >j i he study 01 dark matter halo central densities and mass r ■ . ^ r .1 m • -n ; i.- 1 

5^ . . . . -1 1 1 n At the lamt end 01 the Lh , no signifacant evolution has 

!3 distribution allows us to constrain possible galaxy lornia- , 1 ; , 1 1 , , , 1 1 • i , 11 1 ■ 1 , ■ 

. „ . . been detected, but at the bright end clear brightening was 

■ " ' tion models and large scale structure lormation scenarios „ 1 . ■ ; t . 1 i -n / 1 • -n 

o . r.r.^r. t^, • , A i i r>^^r> louud . biVCu at intermediate redshitts ( z < 1 . 5 ) siguificant 

(Navarro & btemmetz .2000,: Khairul Alam et al. 2002. ,,. ju/^ • 

\, .. ' — ~r 1 . , n 1 • • • evolution was detected at the bright end by Cowie et al. 

Gentile et al. 120041). hoi this kind oi analysis, it is neces- /[nr^nA^ rrn r, , i ,nnnA^ ^ j- j j-i m j- r v 

, ' — ,. . r . ., 1 , , 1 (12004). Toft et al. (2004) studied the LF diagram of an X- 

sary to know both the distribution ot visible and dark mat- 
ter. Unfortunately, the structure and mass distribution 
of stellar populations, and therefore, the discrimination 
between visible and dark matter in galaxies (especially The C-M diagram of a large number of galaxies in the 

in very distant galaxies) is not known precisely enough, redshift range 0.2 < z < 1.1 was studied by Bell et al. 

Thus, it is interesting to analyze the corresponding data l)2()()4|l . and evolutionary trends for both early type and 

at different redshifts. late type galaxies were found. While the early type galax- 

In most cases, the evolution of the visible structure of ies evolve only mildly with redshift, the number of lumi- 

galaxies with redshift has been studied, e.g. the luminos- nous late type galaxies increases remarkably with redshift. 

ity function (LF), C-M diagram, fundamental plane (FP), mi r i ; i i • i / r, ^ 

,1 m 11 XT 1 /rm^\ 1- i i ,i r i , • 1 hc tundamental plane m two clusters (z ~ 0.58 and 

the TuUy-Fischer (TF) diagram, scale-lengths of galactic ^ ^, ^^j^ 

disks. 

The luminosity function in deep fields has been studied 
by Gabasch et al. (2004) in the FORS Deep Field up to 



0.83) was studied by Wuyts et al. H2004|) . They confirm 
an earlier result that the slopes of the FP do not change 
within the studied redshift range. However, there is a gen- 
eral evolutionary shift indicating that galaxies in these 
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To study the evolution of the TF diagram, the maximal 
rotation velocities in rich cluster galaxies at redshifts z — 
0.3 — 0.5 were measured by Ziegler et al. (|2003f) . Their 
analysis did not show any significant difference between 
cluster and field galaxies at the same redshifts, even if 
compared to the local galaxies. At higher redshifts (z — 
0.2 — 1.2, mean z ~ 0.7), the TF diagram was studied on 
the basis of 83 disk galaxies by Conselice et al. (|2004p and 
also no statistically significant evolution with redshift was 
detected. On the other hand, Milvang-Jenses et al. I12003|I 
found that cluster galaxies are brighter by 0.5™ — 1™ in B 
color when compared to field galaxies at the same redshift 
(z = 0.8). For field galaxies at redshifts z — 0.1 — 1.0 
the TF diagram was obtained by Bohm et al. (2004;) and 
brightening of galaxies with redshift depending on galactic 
mass was discovered. 

If, in addition to the surface photometrical data, the 
rotation velocities of galaxies are also available, galactic 
structure can be described in greater detail. In several 
projects, the rotation of distant galaxies has been mea- 
sured. In addition to the works referred in Paper I, Erb et 
al. H2003|l measured rotation curves (RCs) of 6 very dis- 
tant [z > 2) galaxies. Unfortunately, the extent of these 
RCs is rather small (less than 1"), the maximum rotation 
values have often not been reached and the models of these 
galaxies would be uncertain. 

In relation to the study of the TF diagram, Ziegler et 
al. H2003j) measured RCs for 13 galaxies in rich clusters at 
redshifts z = 0.3 — 0.5. Of these galaxies 4 — 5 objects can 
be used for mass modeling (the extent of the measured 
RCs > 2", kinematics more or less regular). The largest 
sample of RCs at intermediate redshifts was measured by 
Jager et al. H20U4|I . where objects in 7 clusters at redshifts 
z = 0.3 — 0.6 were analyzed and several of the derived RCs 
could be used for mass distribution modeling. 

To construct a self-consistent photometrical and dy- 
namical model of a galaxy, both surface photometry and 
rotation curve data are needed. For galaxies at intermedi- 
ate redshifts, the typical scale is 4 — 8 kpc per arcsec and in 
order to determine the parameters of galactic components, 
one should use high-resolution Hubble Space Telescope 
(HST) photometry in addition to the other photometrical 
data. 

In the present study, we continue our modeling of dis- 
tant galaxies (Tamm & Tenjes, [20021 Paper I). To derive 
the surface brightness distribution of these galaxies, we 
used images from the HST archive. RCs were measured by 
Vogt et al. ifTMHI and Rigopoulou et al. General 
properties of these galaxies are given in Table ^ Galactic 
names are from the NED database. In the present work 
we take Hq ~ 65 km s^^Mpc~^ and qo ~ 0.5. 

2. Observations and photometry 

We used the HST WFPC2 images of the galaxies for 
photometry. The images of GSS 104-4024, GSS 064- 
4442 and MDS uemO-043 were retrieved from the HST 
archive, HDFS J223247.66-603335.9 is a galaxy on the 



Hubble Deep Field South, available via the HST web page. 
General data about the observational programs and the 
images are presented in Table 

We used pipeline-reduced, on-the-flight calibrated im- 
ages. The Hubble Deep Field South is available with all 
necessary reduction done, for the other three galaxies we 
conducted exposure combining, cosmic-ray removal and 
background removal with IRAF and STSDAS software 
packages, as described in detail in Tamm & Tenjes 1)20011 
120031 . The resulting images of the galaxies can be seen in 

Figca 

We used the Lucy-Richardson method to deconvolve 
point spreading. We calculated point-spread functions of 
the WFPC2 camera using TinyTim 6.0. The V-hand 
image of GSS 104-4024 could not be deconvolved be- 
cause of too low signal-to-noise and signal-to-background 
ratios. Isophotes were calculated using the ellipse task 
of STSDAS. To construct luminosity profiles from the 
isophotes, we used deconvolved images for the central 
parts of galaxies and original images for the outer radii. 

We calibrated WFPC2 luminosity counts to Johnson 
standard magnitudes using formulae and tables given by 
Holtzman et al l|1995b|l . see also Tamm & Tenjes H2001I 
12003(1 for details. Luminosity was corrected for the cos- 
mological dimming by a factor of (1 -I- z)"*. The resulting 
luminosity profiles of both original and deconvolved im- 
ages in both V and / passbands are presented in Fig. |21 
To present luminosity profiles we prefer a logarithmic scale 
for distances, enabling a better overview of the deconvolu- 
tion results. The scatter of the datapoints of different runs 
of the ellipse task illustrates the uncertainties caused by 
noise and the limited resolution of the images. 

For modeling we derived rest-frame 5-band luminosi- 
ties, applying the k-correction according to van Dokkum 
& Franx (1993|), Kelson et al. (2000). Synthetic spectra 
were taken from the paper by Coleman et al. (1980), the 
relation between Johnson and A_B-magnitudes from Frei 
& Gunn H1994fl . We determined the following transforma- 
tions to be conducted 

B{z = 0.48) = / + 0.585(y - /) -f 0.7 

B{z = 0.58) = / -I- 0.36(V" ~ I) + 0.86 

B{z = 0.81) =1 + 0.95 

B{z = 0.88) = I~ 0.106(1/ - /) + 1.18 

The resulting luminosity profiles in the galactic rest-frame 
B color are shown in Figs |21 and 01 by open circles. These 
profiles were used for further modeling. 

The RCs for the galaxies GSS 104-4024, GSS 064-4414 
and MDS uemO-043, based on Keck spectroscopy, were 
taken from the paper by Vogt et al. H1996|l . The veloc- 
ity profile of HDFS J223247.66-603335.9 was measured by 
Rigopoulou et al. (2002) using VLT spectra. 

We have minimized the difference between the RCs 
from the two sides of each galaxy, shifting the coordinates 
(''0, Vq) of the RC panels (see Persic & Salucci (1995) for 
details). This makes the RC profiles more regular. The re- 
sulting folded RCs are presented with circles in Figs El- El 
We have used different types of circles (open and closed) 
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Name 


RA 


DEC 


Hubble 


2" 


Scale 


Inclin." 


m'} 




M% 




(2000) 


(2000) 


type 




(kpc/") 


(deg) 


(mag) 


(mag) 


(mag) 


GSS 104-4024 


14''17™26.9= 


52° 24' 49. 5" 


Sbc 


0.812 


6.52 


82 


21.88 


0.034 


-20.5 


GSS 064-4442 


14''17™54.0' 


52°29'14.3" 


Sbc 


0.877 


6.62 


55 


22.05 


0.033 


-20.4 


MDS uemO-043 


3''4'"59.1' 


-0°11'47.0" 


Sc 


0.476 


5.51 


46 


19.21 


0.362 


-21.4 


HDFS J223247.66 




-60°33'35.9" 


Sb(c) 


0.581 


5.95 


49 


19.14 


0.119 


-21.2 


-603335.9 





















° Galactic redshifts are from Vogt et al. and (for HDFS J223247.66-603335.9) Rigopoulou et al. iQUO?!! . inclinations 

are from Simard et al. H2002fl and Rigopoulou et al. I|2002^ . 

' Integrated apparent / magnitudes calculated from our photometry. 

Absolute B- magnitudes in the galactic rest-frame are calculated from our models (Sect. 4), considering the absorption in 
the Galaxy Ab according to Schlegel et al. 1199811 . 



Table 2. HST observations used for photometry. 



Name" HST survey total exposure background background 

program time in / (s) level in / level in V 

GSS 104-4024 Groth Strip Survey 4400 5^9 6^9 

GSS 064-4442 Groth Strip Survey 4400 6.15 6.9 

MDS uemO-043 Medium Deep Survey 6700 13.5 23.3 

HDFS J223247.66-603335.9 HDF-South 113 900 0.0 0.0 




Fig. 1. Images of the galaxies (from left to right) GSS 104-4024, GSS 064-4442, MDS uemO-043 and HDFS J223247.66- 
603335.9 through F814W filter after background subtraction. 



to discriminate between the observations in opposite di- 
rections from the center. Error bars are also shown. We 
have tried to select galaxies with kinematics as regular as 
possible. 



3. Model description 

As in Paper I, we limit the main stellar components to 
the bulge and the disk. To construct a dynamical model, 
a dark matter component - the dark halo - must be added 
to visible components. 



To construct a model with two stellar components, the 
surface luminosity distribution is usually approximated by 
a r^/*-bulge and an exponential disk. These surface bright- 
ness distributions can be expressed with the help of the 
general Sersic formula |^ 

/(a) =/(0)cxp[-6„,(aK)i/™], (1) 

(Sersic [T968|l with a parameter m > 2 for the bulge and 
m = 1 for the disk. In this formula a is distance along the 
galactic major axis, Oc is the radius containing half of the 
total luminosity, is a normalizing constant. 




Fig. 2. Surface brightness distribution of the galaxies (from left to right) GSS 104-4024, GSS 064-4442, MDS uemO-043 
and HDFS J223247. 66-603335. 9 in I and V colors in logarithmic distance scale. Both original and deconvolved profiles 
are given. 



If, in addition to the photometrical data, kinematic 
data are also used, the corresponding dynamical model 
must be consistent with the photometry, i.e. the same den- 
sity distribution law must be used for RC modeling (and 
for the velocity dispersion curve, if possible). In the case 
of an infinitely thin disk, the circular velocity for an expo- 
nential surface density distribution can be expressed via 
modified Bessel functions fFreeman ll970(l . In addition, for 
spherical systems an expression for circular velocity with 
an integer Sersic index m (Mazure & Capelato l2UU2|l can 
be derived. For a non-integer index and ellipsoidal surface 
density distribution a consistent solution for RC calcula- 
tions is not known. 

The density distribution parameters are often deter- 
mined by the least squares method. In this case the pa- 
rameter m is not necessarily an integer number. In the 
present paper, the density distribution parameters are de- 
termined by the least squares method and they can have 
any value. Moreover, we have found that at least for our 
present sample of high redshift galaxies the best fit for the 
disk luminosity distribution is achieved with the Sersic in- 
dex m < 1. 

For the reasons given above we decided in addition to 
models with a r^/** bulge and exponential disk to construct 
models for the galaxies starting from a spatial density dis- 
tribution law, which allows an easier least-square fitting 
simultaneously for light distribution and RC. 



In such models, the visible part of a galaxy is given as 
a superposition of a bulge and a disk. The spatial density 
distribution of each visible component is approximated by 
an inhomogeneous ellipsoid of rotational symmetry with 
the constant axial ratio q and the density distribution law 

p(a)=p(0)exp[-(a/(fcao))i/^], (2) 

where p(0) = hM/{47TqaQ) is the central density and 
M is the component mass; a = y'R^+z^Jq^, where R 
and z are two cylindrical coordinates, oq is the harmonic 
mean radius which characterizes rather well the real ex- 
tent of a component, independently of the parameter N. 
Coefficients h and k are normalizing parameters, depend- 
ing on N, which allows the density behavior to vary with 
a. The definition of the normalizing parameters h and k 
and their calculation is described in Tenjes et al. 1)1994(1 . 
Equation allows a sufficiently precise numerical inte- 
gration and has a minimum number of free parameters. 

The dark matter (DM) distribution is represented by 
a spherical isothermal law 

p(0)([l + (t)T^-[l + (£)T^) «<«° (3) 
lO a>a°. 

Here a° is the outer cutoff radius of the isothermal sphere, 
flc = kao- 

The density distributions for the bulge and the disk 
were projected along the line of sight, divided by their 
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Fig. 3. Surface brightness distribution of the galaxies GSS 104-4024, GSS 064-4442, MDS uemO-043 and HDFS 
J223247. 66-603335. 9 in rest-frame B color (open circles). Dashed lines - surface brightness distribution of best fit 
r^/* bulge and exponential disk, solid line - total surface brightness distribution. 



mass-to-luminosity ratios / and their sum gives us the 
surface brightness distribution of the model 

«^>-|:&/;<^' 

where A is the major semiaxis of the equidensity ellipse of 
the projected light distribution and Qi are their apparent 
axial ratios = cos^ 7 + 9^ sin^ 7. The angle between the 
plane of a galaxy and the plane of the sky is denoted by 
7. The summation index i designates two visible compo- 
nents, the bulge and the disk. 

The masses of the components were determined from 
the rotation law 

„ pj(a)a^da , , 

^nR)-^^,.Gi ^J^^^i, . (5) 

V\R) = Y^vl{R), (6) 



where G is the gravitational constant, e = \J\ ~ (p' is ec- 
centricity, and R is the distance in the equatorial plane of 
the galaxy. Now the summation is over all three compo- 
nents. 

We have kept the axial ratio of the components fixed, 
taking q — 0.7 for the bulge and q = 0.1 for the disk, ac- 
cording to analogy with the near-by galaxies. The model 
parameters oq, Lb, M and N for the bulge and the disk 
were determined by a subsequent least-squares approxi- 
mation process. First, we made a crude estimation of the 
population parameters. The purpose of this step is to avoid 
obviously non-physical parameters - relations (0} and © 
are nonlinear and fitting of the model to the observations 
is not a straightforward procedure. Next, a mathemati- 
cally correct solution was found for each galaxy. Details 
of the least squares approximation and the general model- 
ing procedure were described by Einasto & Hand H1989fl . 
Tenjes et al. (fl^MlIT^ . 

The RCs of the modeled galaxies have a rather small 
extent (see Figs.EHHl. We do not have any additional in- 
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Fig. 4. Surface brightness distribution of the galaxies GSS 104-4024, GSS 064-4442, MDS uemO-043 and HDFS 
J223247. 66-603335. 9 in rest-frame B color (open circles). Dashed lines - surface brightness distribution of best fit 
model components using spatial density distribution Solid line - total surface brightness. 



formation about the dynamics of these galaxies at large 
galactocentric radii. For this reason, the outer cutoff ra- 
dius of the DM component remains undetermined at 
present. As in Paper I, we fixed = Sag on the basis of 
the structure of nearby galaxies. This value influences the 
behavior of the RC only in the extreme outer regions and 
does not influence our results. In addition, by analogy with 
nearby galaxies we assumed that the RCs of these galax- 
ies remain nearly flat at least up to 30 kpc ~ 1.5 — 2Ropt 
(Persic et al lTtMl Sofue & Rubin IMnil . 

As we have no information about velocity dispersions 
in the central regions of our galaxies, we can calculate only 
the circular velocities (Eqs. [SlEl). The difference between 
the circular and rotational velocities is known in galactic 
dynamics as the asymmetric drift and is mainly a function 
of velocity dispersions. Typical emission-line dispersions in 
disk galaxies at intermediate redshifts are 30-100 km/s or 
even more (Im et al. lMTTl Erb et al. 11003) • Thus, within 
the central 0.5" — 0.7", modeled velocities must be higher 
than the observed rotational velocities (see also Sect. 4). 



4. Results 

We have fitted the two-component models of the galaxies 
to the B color profiles. Profiles fitted by a r^/^-bulge and 
an exponential disk can be seen in Fig. O Profiles fitted 
on the basis of the spatial density law 10) are presented 
in Fig. 01 In both cases dotted lines correspond to the in- 
dividual components - the bulge and the disk. The solid 
line is the sum of the two components. It is seen that a 
non-exponential disk gives a better agreement with ob- 
servations. The mean deviations of the models from the 
observations in Fig. Hare only - /im°d<=i) = 0.07'", 

0.04", 0.14" and 0.13™, respectively. 

Fitting of the modeled RCs to the observations is much 
more complex. First, the resolution of the present-day 
spectroscopic observations is insufficient to construct high- 
precision velocity profiles for galaxies at such high red- 
shifts. In addition to the uncertainties and the low extent 
of the observations, the velocities in the central regions 
remain severely underestimated due to the lack of reso- 
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lution and therefore the bulge mass remains practically 
unknown. 

It is possible to fit only the share of the disk or the 
dark matter component to the observed rotation and still 
achieve satisfactory concordance. In order to keep the re- 
sults more realistic we have, however, constructed three- 
component models for all four galaxies with a bulge, a 
disk and a dark matter component. For the reasons given 
above, the mass of the bulge component should not be 
taken too literally; we have practically ignored the ob- 
served velocities at low radii. As in most cases no other 
hint of the mass of the bulge was given, we took the 
mass-to-light ratio of the bulge equal to that of the disk. 
Assuming the bulge of a galaxy to be populated by stars 
from an older generation than the disk, we should thereby 
get the approximate lower limit of the bulge mass. 

We have scaled the original observed RCs to match 
"l^nax given by the authors, thereby compensating for the 
effects of misalignment of the slit, "beam smearing" and 
seeing at the outer radii, which we have used to fit our 
modeled RC. 

The parameters of the DM component cannot be de- 
termined uniquely by our method. We could, for example, 
increase its mass, if we increased r simultaneously. A pa- 
rameter that can be determined reasonably well, however, 
is the central density of the dark matter p(0) (see Sect. 5). 

4.1. Notes about individual galaxies 

4.1.1. GSS 104-4024 

This Sbc galaxy at z = 0.8116 is seen nearly edge-on (in- 
clination 82°); a possible dust lane along its semi-major 
axes can be detected on the HST image (Fig.^, obscuring 
the central part of the bulge. The usual isophote fitting 
would underestimate the luminosity of the central regions 
(r < 0.3"); the true luminosity remains unknown. A possi- 
ble correction suggested by our modeling software is shown 
in Fig. 0](the correction contributes only 0.1™ to the total 
luminosity) . 

We calculated the total apparent V-band luminosity 
of the galaxy to be mj — 21.88 ± 0.06, which matches 
well with the estimate 22.05 by Simard et al. (2002). The 
mean apparent color index {V ~ I) — 2.1 ± 0.4. Taking 
into account the absorption in the Galaxy, we found the 
absolute B-band luminosity to be Mb = —20.5. However, 
considering also the internal absorption, the galactic lu- 
minosity would be higher and the resulting M/ L ratio of 
visible matter would be smaller. For the exponential bulge 
we find = 11 kpc, which coincides with the value 11.54 
given by Simard et al. plH^ . 

4.1.2. GSS 064-4442 

This is an Sbc galaxy aX z — 0.877, the most distant one 
discussed in this paper. We calculated the total apparent 
Johnson m/ = 22.05 ± 0.14, which matches well with the 
estimated value 21.97 by Simard et al. ((2002,1 . The mean 



apparent color index {V — I) — 1.5 ± 0.2, absolute B- 
band magnitude Mb = —20.4. We obtained = 2.2 kpc 
for the exponential disk, which somewhat differs from the 
estimate 3.76 by Simard et al. H2()()2|l . Taking into account 
the faintness and low resolution of the image of the galaxy 
this deviation is not surprising. 

4.1.3. MDS uemO-043 

This well-exposed Sc galaxy at z = 0.476 has an obvious 
spiral structure. We found the total Johnson /-band lu- 
minosity to be 19.21 ± 0.04, which agrees well with the 
result of Vogt et al. 1)199611 mj = 19.31. The mean ap- 
parent color index {V — I) = 1.5±0.2. With correction 
for the absorption in the Galaxy applied, MDS uemO-043 
turns to be a bright galaxy: Mb = —21.4. We estimated 
rd — 4.0 kpc for the exponential disk, which is reasonably 
close to rd = 5.24 by Vogt et al. H1996|l . However, in Vogt 
et al. H1996|l deconvolution of the surface brightness pro- 
file has not been done, thus a slight difference in disk scale 
is normal. 

4.1.4. HDFS J223247.66-603335.9 

This bright and massive Sb galaxy is very well exposed 
on the HDF-South image and was studied in a paper by 
Rigopoulou et al. (2002). The excellent quality of the HDF 
images makes it possible to analyze the photometry of the 
galaxy with high precision. We calculated the total appar- 
ent Johnson / luminosity to be 19. 14 ±0.02. Transforming 
this result to the AB system (adding 0.35 mag, Frei & 
Gunn,E^, we get a perfect match with the estimate of 
19.5 by da Costa et al. H1998|l on the basis of ground-based 
observations with the VLT. The mean apparent color in- 
dex {V — I) — 2.0±0.2. The absolute B-band luminosity is 
AIb = —21.2. We obtained rd = 3.3 kpc for the exponen- 
tial disk, which unfortunately cannot be compared to the 
value by Rigopoulou et al. |2002f) due to different surface 
brightness scaling. 

The velocity profiles are presented in Fig.|Sl Calculated 
from our model, the total mass within 28 kpc is 13 • 
IO-'^^Mq, which agrees well with the mass 10 • IO^^Mq 
within 20 kpc estimated by Rigopoulou et al. I)2()()2|l . 

4.2. Calculated model parameters 

In Table 121 calculated from the models, the integrated pa- 
rameters of a Sersic surface density distribution are given: 
the scale lenghts of the bulge ri, and the disk rd and the 
corresponding central luminocities in rest-frame -B-band 
(/i-^)buigc, (M'B)disk- The corresponding profiles can be seen 
in Fig. ill 

The parameters of the final models on the basis of the 
spatial density distribution (0) (the harmonic mean radius 
flo, the total mass of the population M, the structural pa- 
rameter N, the dimcnsionless normalizing constants h and 
k, B-luminosities and the corresponding mass-to-light ra- 
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Table 3. Derived parameters for a Sersic surface density 

distribution. q 5 



r (kpc) 
10 



15 



Name" 


^6 






(Ms)disk 




(kpc) 


(kpc) 


(mag) 


(mag) 


GSS 104-4024 


2.0 


11.0 


14.0 


21.2 


GSS 064-4442 


0.3 


2.2 


14.5 


19.8 


MDS uemO-043 


0.23 


4.0 


13.0 


20.4 


HDFS J223247.66 


0.3 


3.3 


12.7 


20.2 


-603335.9 











Table 5. Derived integrated parameters from models. 



Name 


IO^^Mq 


Lb 

10i»L« 


M/Lb 


D/B 


GSS 104-4024 


11. 


2.53 


4.3 


4.0 


GSS 064-4442 


2.1 


2.34 


0.9 


6.5 


MDS uemO-043 


7.6 


5.43 


1.4 


12.2 


HDFS J223247.66 


34. 


4.58 


7.4 


5.1 


-603335.9 
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Fig. 5. Rotation curve of the galaxy GSS 104-4024. Open 
and gray circles - observed velocities by Vogt et al. (1996) 
at the opposite sides of the galaxy, thick solid line - rota- 
tion curve from the mass model, dashed lines - rotation 
velocities due to visible components, thin solid line - ro- 
tation velocities due to the dark matter component. 



tics) are given in Table 0] The masses of the bulge and 
DM are given in parentheses to indicate that on the basis 
of the present observational data it is not possible to de- 
termine these values uniquely. The final model profiles are 
denoted by thick solid lines in Figs. Models of indi- 
vidual visible components are presented by dashed lines, 
the model of the dark matter component by a thin solid 
line. 

In Table |5l the total mass of the visible matter Mvis, 
the total intrinsic luminosity Lg, the mass-to-light ratio 
of the visible matter M/Lb and the disk-to-bulge ratio 
D/B are given. 
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Fig. 6. Rotation curve of the galaxy GSS 064-4442. 
Designations are the same as in Fig. O 
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Fig. 7. Rotation curve of the galaxy MDS uemO-043. 
Designations are the same as in Fig. |5| 
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Fig. 8. Rotation curve of the galaxy HDFS J223247.66- 
603335.9. Designations are the same as in Fig. [S] 

5. Discussion 

We have used the isothermal approximation for the dark 
matter density distribution in our models. However, sev- 
eral other dark matter distributions have been suggested 
and used by different authors. The main debate has been 
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Table 4. Model parameters for spatial density distribution. 
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Name 


ao 


M 


M/Lb 


Lb 


iV 


h 


k 




(kpc) 


(IO^OMq) 


(Mq/Lq) 










GSS 104-4024 
















bulge 


1.5 


(2.2) 


4.3 


0.50 


0.91 


3.33 


0.594 


disk 


13.2 


8.7 


4.3 


2.03 


0.32 


1.18 


1.37 


dark matter 


80. 


(550.) 








14.8 


0.151 


GSS 064-4442 
















bulge 


0.9 


(0.3) 


0.90 


0.31 


0.51 


1.59 


1.12 


disk 


5.0 


1.8 


0.90 


2.03 


0.34 


1.23 


1.34 


dark matter 


50. 


(230) 








14.8 


0.151 


MDS uemO-043 
















bulge 


0.5 


(0.6) 


1.41 


0.41 


1.11 


4.96 


0.405 


disk 


7.2 


7.1 


1.41 


5.02 


0.73 


2.38 


0.805 


dark matter 


80. 


(340.) 








14.8 


0.151 


HDFS J223247.66 
















-603335.9 
















bulge 


0.7 


(5.6) 


7.4 


0.75 


0.80 


2.70 


0.719 


disk 


7.0 


28.3 


7.4 


3.83 


0.67 


2.12 


0.888 


dark matter 


80. 


(1500.) 








14.8 


0.151 



between "cuspy" distributions deduced from N-body sim- 
ulations, with the density increasing infinitely towards the 
center of the halo, and "cored" distributions usually de- 
rived from the observed RCs of late type and low surface 
brightness galaxies with a constant density distribution in 
the central part of the halo. 

According to their numerical simulations, Navarro, 
Frenk & White fl996, NFW) proposed the following for- 
mula for the density distribution 



p(r) 



{r/rs)[l + {r/rs)Y 



(7) 



where ps and Tg are the characteristic density and radius, 
respectively. Moore et al. H1999|l proposed a rather similar, 
but still "cuspier" distribution: 



p(r) 



PsM 



(r/r,)i-5[l + (r/r,) 



1.51 



(8) 



Burkert fl995l) studied low surface-brightness galaxies 
and approximated the dark component with a distribution 
similar to isothermal 



p{r) 



(r + ro)(r2 + r^) 



(9) 



We have constructed a "best estimate" curve for the 
velocity distribution of the galaxy MDS uemO-043, using 
the observed rotation for the inner radii and keeping the 
rotation velocities beyond them constant. Thereafter we 
subtract the share of the modeled visual components from 
the resulting curve, revealing an estimate for the share of 
the DM component V^y^ — V^bs~K?is- These velocities are 
shown as open circles in Fig. [Altogether with the estimated 
errors and they correspond to the supposed rotation con- 
tribution of DM. Finally, we test how well different density 
distributions fit with the derived DM contribution. The 
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Fig. 9. Comparison of the dark matter contribution in 
kinematics of the galaxy MDS uemO-043 (open circles with 
estimated error bars) with the different dark matter den- 
sity distribution models. Two vertical arrows indicate out- 
ermost observed points of the surface brightness profile 
and the rotation curve. 



large uncertainties at small radii are caused by the effects 
of limited spectral and spatial resolution. 

Except for the "cuspiest" one by Moore et al. H1999|l . 
the curves of all other distributions fit within the error 
bars (Fig. Burkert 's model is nearly the same as our 
isothermal law in the central part, while the best fit of 
NFW stays slightly above these two at small radii and 
below them at higher radii. Due to the limited spatial 
resolution of the RCs of high redshift galaxies, it is not 
possible to prefer the NFW, Burkert or isothermal distri- 
butions, nor to study the most recent profiles by Navarro 
et al. H2004(l and Fukushige et al. I|2004|l . We intend do 
this analysis on the basis of detailed models of nearby 
well-studied galaxies. 
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As the behavior of the density distributions (7) and 
(8) at the center is limited by the resolution of N-body 
calculations and the behavior of density distributions (3) 
and (9) is quite similar (Fig.|HJ), below we discuss only the 
isothermal DM halo. 

Due to the limited extent of RCs it is not possible 
to determine separately the radius and the mass of the 
DM distribution. However, the combination of these two 
parameters, the central density of the dark matter p(0) = 
/iM/(47rao), characterizing the central rotation velocity 
gradient caused by the dark halo, is nearly independent 
of the choice of the dark halo outer cutoff radius and can 
be determined. Calculated from our models, the central 
densities of DM components are p(0) = 0.013, 0.022, 0.008 
and 0.035 in units Mq/pc^. These values are calculated 
for maximum disk models and are thus the lower limits. 
Adding the two galaxies at redshifts 0.9 and 0.99 from 
Paper I we get the range of the DM central density of 
the four galaxies at the mean readshift (z) ~ 0.9 is p(0) — 
(0.012-0.028) Mq/pc^. This value can be compared with 
the mean central DM density of nine local galaxies p{0) = 
(1 - 4) • 10-24g/cm3 = (0.015 - 0.050) Mg/pc^ derived 
by Borriello & Salucci H2001|l . Although our sample of 
galaxies is small, there seems to be no significant evolution 
of DM central density with redshift. 

We also calculated the amount of DM within a half- 
light radius. For galaxies at redshifts (z) ~ 0.9 the mean 
dark matter contribution was 40 ± 15 per cent. Thus in 
our case, the discrepancy between the observed galaxies 
and CDM models is not as high as was indicated in the 
analysis by Bell et al. H2003|l . 

According to our models, the mass-to-light ratios of 
visible matter (within the maximum disk assumption) are 
M/Lb = 4.3, 0.9, 1.4 and 7.4. Adding the two galaxies 
at redshifts 0.9 and 0.99 from Paper I to these galaxies, 
we find that the mean M/Lb of the four galaxies at the 
mean redshift {z) ~ 0.9 is {M/Lb) = 2.5. Another way to 
determine stellar masses and mass-to-light ratios is from 
broadband optical-NIR photometry and chemical evolu- 
tion models (Drory et al. fi200'l)) . Berta et al. H2004|l ). For 
galaxies at redshifts z ~ 0.9 Drory et al. H2004|l derived 
{M/Lb) ~ 1-3, being somewhat smaller than the value 
calculated by us. For the highest mass galaxies at redshift 
z ~ 0.5, they derived {M/Lb) = 2.6, also smaller than 
our value of 7.4 for the galaxy HDFS J223247.66-603335.9. 
However, we think that within the precision of the mod- 
els used these results agree reasonably. Stellar masses of 
starburst galaxies at z = 0.5 are found to be even as large 
as 50 • IO^^Mq (Drory et al. (2004), Berta et al. (2004,) ). 

The derived M/Lb = 2.5 at z ^ 0.9 can be compared 
to the mean stellar M/Lb calculated by Giraud ('1998| 
for local Sb galaxies M/Lb = 2.7. No significant evolu- 
tion with redshift is detected - a result which was found 
earlier in the studies of the TF diagram at intermediate 
redshifts (Vogt et al. 119961 Ziegler et al. i2003i Conselice et 
al. H2004|l . Bohm et al. H2004L for higher mass galaxies). An 
explanation for these small changes can be found in chem- 
ical evolution models. Evolution calculations by Bicker et 



al. H2004L Fig. 7) indicate that the k-correction and evo- 
lution corrections in the i?-band nearly compensate each 
other in the case of Sb galaxies at redshifts z ~ 0.5 — 1. 

As mentioned in Paper I, the disk concentration pa- 
rameter N seems to have a trend towards lower values at 
higher redshifts. Even visual inspection of the luminosity 
profiles (Figs [5]- 0J reveals that the surface brightness 
of the outer parts of the galaxies decreases faster than ex- 
ponentially. Surface brightness distribution of disks with 
the concentration parameter A'^ < 1 is similar to the sur- 
face brightness distribution of truncated disks described 
by Pohlen et al. ("2000, 2004). According to recent SPH 
simulations, truncated disks seem to be a natural phenom- 
ena within ACDM models (Governato et al. 120041 Fig. 2). 
However, in our case the truncation takes place generally 
slightly earlier than at 3-5 disk scale lengths. Truncation 
of a sample of redshift z ~ 1 disk galaxies has recently 
been studied by Perez (2004) . 

It is possible to construct mass distribution models for 
galaxies without assuming the existence of a DM compo- 
nent. The corresponding models are given in Fig 1101 Now 
the calculated RCs start to decrease significantly near or 
right after the last observed point of rotation. The mass- 
to-light ratios for the disks are in this case M/Lb = 12., 
1.8, 1.9 and 9.7, respectively. Thus, within the presently 
available observations, the models without any dark halo 
fit well in most cases with the observations and it is not 
possible to determine the existence of a DM component in 
these four galaxies on the basis of their kinematics. The 
mean mass-to-light ratio for z ^ 0.9 galaxies (adding the 
two galaxies from Paper I) is now {M/Lb) = 8.5. 

Figure EH shows that RCs can be well fitted without 
DM. However, we think that there are enough independent 
arguments in support of the existence of DM and as a 
result, we think that realistic models must include a DM 
component. 
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